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ABSTRACT 

We present high precision transit observations of the exoplanet WASP-21b, obtained with the 
RISE instrument mounted on 2.0m Liverpool Telescope. A transit model is fitted, coupled 
with an MCMC routine to derive accurate system parameters. The two new high precision 
transits allow to estimate the stellar density directly from the light curve. Our analysis sug- 
gests that WASP-21 is evolving off the main sequence which led to a previous overestimation 
of the stellar density. Using isochrone interpolation, we find a stellar mass of 0.86 ± 0.04Mq 
which is significantly lower than previously reported (l.Ol ± O.OSM© ). Consequently, we 
find a lower planetary mass of 0.27 ± O.OlMjup. A lower inclination (87.4 ± 0.3 degrees) 
is also found for the system than previously reported, resulting in a slightly larger stellar 
(i?, = 1.10 ± 0.03 R0) and planetary radius = 1.14 ± 0.04 Rjup). The planet radius sug- 
gests a hydrogen/helium composition with no core which strengthens the correlation between 
planetary density and host star metallicity. A new ephemeris is determined for the system, i.e., 
0"= 2455084.51974 ± 0.00020 (HJD) and P = 4.3225060 ± 0.0000031 days. We found no 
transit timing variations in WASP-21b. 
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1 INTRODUCTION 

Transiting planet systems are valuable because their geometry en- 
ables us to estimate accurate planetary properties. Time-series pho- 
tometry during the transit allows us to derive the orbital inclina- 
tion and the relative radii of the host star and planet. Combining 
this with radial velocity variations and stellar parameters, allows us 
to derive the absolute mass of the planet. Hence, the bulk density 
of the planet can be estimated with good accuracy, giving us in- 
sight into its composition (Guillot 2005; Fortney et al. 20071), thus 
placing constraints on planetary structure and formation models. 
Given the remarkable diversity in the structure of large planets, it is 
important to obtain planetary parameters which are as accurate as 
possible. However, obtaining high signal-to-noise transit observa- 
tions is difficult and consequently even some of the brightest stars 
with planets are lacking good quality light curves and, hence, have 
poorly determined planetary parameters. 

The RISE (Rapid Imager to Search for Exoplanets) instru- 
ment, mou nted on the 2.0m Liverpool telescope dSteele et alj2008l ; 
iGibson eT a l. 2008) was designed for exoplanet transit observa- 
tions. Its main scientific driver was the detection of transit-timing 
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variations and hence the search for low-mass companions to "hot 
lupiters". RISE has a rapid readout frame transfer CCD and in 2x2 
binned mode has a readout time of less than 1 s. This implies that 
for exposures longer than 1 s, dead time is negligible substantially 
increasing the time on target. However, most exoplanet host stars 
are relatively bright and saturate the CCD for 1 second exposure. 
To avoid dead time losses, RISE observations are always defo- 
cussed (e.g. Gibson et al. 2008; loshi et al. 2009). Defocussed pho- 
tometry observations have also the advantage of spreading the PSF 
over a larger number of pixels, thereby decreasing flat-fielding er- 
rors. RISE is therefore ideal for obtaining high quality transit light 
curves for exoplanets. 

WASP-21b is a Saturn-mass planet with Mp = 0.30 ± 0.01 
Mjup in a 4.3 day circular orbit (Bou chy et al„ 2010). Its host star 
is a G3V type with M, = 1.01 ± 0.03 Mq ,T^ff = 5800 ± 100 K 
and a low metallicity, [M/H ] =—0.4 ± 0. 1. It was discovered by the 
SuperWASP-North s urvey jPoUacco et al.|[200^ in its 2008-2009 
observing campaign. i BouchvetaUfeOlOh argue that WASP-21 is 
a member of the galactic thick disc because of its low metal abun- 
dances, velocity relative to the Sun and age ~ 12Gyr, which are 
similar to the thick disc population. WAS P-21b is among the low- 
est density planets, pp = 0.24 ± 0.05pj teouchv et alj2010h , and 
has one of the lowest metallicity host stars. Therefore, its prop- 
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erties are particularly importan t for irradiation mo dels. The cur- 
rent parameters of the system teouchv et aTll2010l) are based on 
the SuperWASP discovery photometry and a partial transit light 
curve taken with RISE. However, the lack of a high precision com- 
plete transit light curve required the assumption of the main se- 
quence mass-radius relation which tends to bias the estimate of the 
inclination. Furthermore, the age derived for WASP-21 is longer 
than the main-sequence life time of a I.OIM© star. This suggests 
that WASP-21 could be evolved which would invalidate the main- 
sequence assumption and bias the parameters of the system. To test 
the main-sequence assumption we obtained further observations of 
WASP-21. 

In this paper, we present transit observations of WASP-21b 
with RISE including a full transit light curve. Our high precision 
light curves allow us derive the planetary and stellar radii without 
assuming the main-sequence mass-radius relation for the host star. 
We describe our observations in Section 2. In Section 3, we discuss 
our transit model and present the updated parameters of the sys- 
tem in Section 4. Finally, we discuss and summarise our results in 
Section 5. 



2 OBSERVATIONS 

WASP-21b was observed with RISE ^Steele et al.l2008h mounted at 
the auxiliary Cassegrain focus of the robotic 2.0m Liverpool Tele- 
scope on La Palma, Canary Islands. This is a focal reducer system 
utilizing a frame transfer e2v CCD sensor. The detector has a pixel 
scale of 0.54 arcsec/pixel that results in a 9.4 x 9.4 arcmin field of 
view. RISE has a wideband filter covering ~ 500-700 nm which 
corresponds approximately to V-l-R. The instrument has no moving 
parts. 

The Liverpool Telescope has a library of flat fields which are 
taken manually every couple of months. RISE fiats are taken during 
twilight at different rotator angles so that there is a uniform illumi- 
nation of the CCD. The exposure times of the images are automat- 
ically adjusted so that the peak counts in the individual flats are be- 
low the non-linearity limit of the CCD at 45000 counts. Typically, 
the individual flats have between 20000 and 40000 counts. Due to 
the fast readout, we can obtain approximately 200 flat frames in a 
run, these are combined to create a master flat. For each observa- 
tion run we use the master flat that is closest in time, although we 
note that these are very stable. 

On 2009-09-09 we obtained a full transit of WASP-21b. A 
total of 6581 exposures in the 2x2 binning mode with an exposure 
time of 2.7 seconds were taken. The telescope was defocussed by 
-1.2mm which resulted in a FWHM of ~ 11 ". For defocussed 
photometry, the star profiles are not Gaussian. However, we found 
that, in our case, a Gaussian provided a good fit to the wings of 
the star profile, and could be use as a rough estimate of the profile 
width. Therefore, we estimated the FHWM in the usual way by 
cross-correlating a Gaussian profile with that of the star. 

A second full transit observation of WASP-21b was attempted 
on 2010-11-24. In this case, deteriorating weather terminated the 
observations shortly after the mid-transit, by which time, 4008 inte- 
grations had been obtained. During these observations, the FWHM 
was ~ 12.5 ". 

Both data- sets were reduced using the ULTRACAM pipeline 
jDhillon et al]|2007h which is optimized for time-series photome- 
try. Initially, we bias subtracted the data while we investigated sys- 
tematic effects that were introduced by the flat fielding process. 
We performed differential photometry relative to five comparison 
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Figure 1. Phase-folded RISE light curves for WASP-21. From top to bot- 
tom in chronological order; 2008 October 07, 2009 September 09 and 2010 
November 24. We superimpose the best-fit transit model and also show the 
residuals for each light curve at the bottom of the figure. The data are binned 
into 30 second periods, and bins displaced vertically for claiity. The indi- 
vidual RISE fight curves pfotted here are avaiiabfe in efectronic form at 
CDS. 

stars in the field, confirmed to be non variable, and we sampled dif- 
ferent aperture radii and chose the aperture radius that minimised 
the noise. For the first night, we used a 22 pixel aperture radius 
(~ 12"), and for the second transit, a 32 pixel aperture radius 
(~ 17"). The photometric errors include the shot noise, readout 
and background noises. 

We also included in our analy sis, the previously published 
egress of WASP-21 taken with RISE dSouchv et al.l20IOl) . For con- 
sistency, we re-reduced the original data using the same method as 
for the other two observations. On 2008-10-07, 2220 exposures of 
5 sec duration were taken. We estimated a FWHM of '--^ 2.7 ", 
therefore, the level of defocussing was lower than in our observa- 
tions. The best aperture radius was found to be 15 pixel (~ 8"). 
Our results agree well with the previous published light curve. 

The final high precision photometric light curves are shown in 
Figure[T]along with the best-fit model described in Section [331 We 
overplot the model residuals and the estimated uncertainties which 
are discussed in Section [T2l 

2.1 Optimum exposure time for RISE 

As mentioned above, defocusing is commonly used in exoplanet 
transit observations. Southworth et al. ( 2009) calculated the opti- 
mum exposure time for the DFOSC imager mounted on the 1.54m 
Danish Telescope. We follow the same procedure and apply it to 
RISE mounted on the Liverpool Telescope and hence, we account 
for readout noise, photon, background and scintillation noise. Sim- 
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ilar to I South worth et ( l2009l) . we do not include flat-fielding 
noise, assuming that the profile position is stable. 

The key difference is that RISE is a frame transfer CCD whose 
dead time is the frame transfer time, 35 milliseconds for observa- 
tions longer than 1 second. For the brightest comparison star in our 
field, (V ~ 9), we found the optimum exposure times with RISE 
are approximately 2.7, 7.8, 10.8 seconds during bright, gray and 
dark time, respectively. 

We iterate that the improve ment in signal-to-noise f or defo- 
cussed observations, reported by ISouthworth et al.l j2009b is only 
due to deadtime losses; hence, the defocussing needed is propor- 
tional to the CCD readout time. If the deadtime was zero the best 
theoretical signal-to-noise would always be for focused observa- 
tions, mainly due to the increase in background noise for wider 
profiles. 

Moreover, in our case, the improvement on signal-to-noise be- 
tween 1 second and 10.8 seconds exposure times is quite small 
on the order of lOppm per 30 sec bin. As we will see below, the 
strongest reason for defocussing is to minimise systematic noise 
which, due to its nature, is not accounted for in the calculation and 
can substantially increase the noise in a transit light curve. Figure|2] 
shows systematic noise variations larger than 400 ppm. 
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Figure 2. Light curve of the brightest comparison star for the 2009 Septem- 
ber 09 observation relative to the ensemble of comparison stars used in 
WASP21b final light curve. It shows systematic noise with an amplitude of 
400 ppm. This compaiison was not used in the final light curve of WASP- 
21. We also overplot the photomefiic errors. 



instrument field of view. We also improved the telescope guiding 
system's stability. This led to an improvement in the precision of 
the light curves which is evident in the latest light curve of WASP- 
21 taken after the upgrades (see Fig. 1). In the November 2010 
observations the variation in position is less than 2 pixels in the x 
direction and 4 pixels in the y. 



3 DATA ANALYSIS 
3.1 Systematic noise 

Exoplanet transit observations are often dominated by systematic 
noise. Therefore, to improve the precision of the light curves it 
is important to determine and minimise this noise source. For the 
2009 September 09 observations, the brightest comparison star (cl) 
on the field was affected by systematic noise. This can clearly be 
seen in Figure (2] where we show the flux of cl relative to the en- 
semble of comparison stars used in the final 2009 WASP-21 light 
curve. This shows a variation of 400 ppm. We found that this sys- 
tematic noise was correlated with the star position in the CCD, 
which during the transit observation varied by 10 pixels in the x 
direction and 8 in the y. Given that we used an aperture radius of 
22 pixels, this implies that only half of the pixels used to perform 
aperture photometry were common for the duration of the obser- 
vation. Hence, we concluded that the systematic noise was due to 
variations in the pixel-to-pixel sensitivity which were not corrected 
by flat fielding. In fact, the systematic noise is slightly higher if we 
flat field the data. Our master flat is a combination of 150 frames, 
each with a mean of 35000 counts. The uncertainty in this flat is 
0.5 millimags per pixel which is smaller than the photometric error 
(~ 4.4 milimags per unbinned point) and the observed systematic 
noise. After careful analysis of the data, we found that the cl com- 
parison star crossed a reflection feature in the CCD that is rotator 
dependent (LT is on an alt-azimuth mount) and thus was not cor- 
rected by flat fielding. This experience demonstrates the importance 
of good guiding in decreasing the sources of systematic noise. If the 
observations were performed in focus and assuming the seeing was 
I arcsec, the FHWM would have been ~ 2 pixels. Using an aper- 
ture radius of 1.5 x FWHM — 3 pixels, it would have implied that 
there were no common pixels during the observations. Therefore, 
we infer, if the observations were focused, the amount of system- 
atic noise would have doubled. Note that the defocussing does not 
affect the guiding since the guide camera is always kept in focus. 

After this incident the RISE instrument was upgraded. The 
source of the reflected feature was identified and removed from the 



3.2 Photometric errors 

An accurate estimate of the photometric errors is important to ob- 
tain reliable system parameters. Our first estimate of errors for each 
light curve includes only the shot noise, readout and background 
noise, which underestimates the true errors. To obtain a more re- 
liable estimate we begin by scaling the errors of each light curve 
so that the reduced of the best fitting model is 1.0. This re- 
sulted in the multiplication of the errors by 1.97, 1.22 and 1.44, 
for the 2008, 2009 and 2010 light curves, respectively. We then 
calculated the tim e-correlated noise following the procedure from 
iGiUon et al. I (12009). Using the residuals of the best fit model, we es- 
timated the amplitude of the red noise, Ur to be 150 ppm, 250 ppm 
and 150 ppm, for the 2008, 2009 and 2010 light curves, respec- 
tively. These were added in quadrature to the rescaled photomet- 
ric errors and were used in the final Markov Cha in Monte Carlo 
(MCMC) chains. However. ICarter & Wind ^2009 !) found that this 
"time-averaging" method of estimating the correlated noise can still 
underestimate the uncertainties by 15-30 per cent. 

3.3 Determination of system parameters 

To determine the planetary and orbital parameters, we fitted the 
thr ee RISE light cu rves of WASP-21b simultaneously. We used 
the lMandel & Agoli j2002l) transit model parametrised by the nor- 
malised separation of the planet, a/7?, , ratio of planet radius to star 
radius, Rp/R,, orbital inclination, i, and the transit epoch, To, of 
each light curve. Our model was originally developed to measure 
transi t timing variations of exoplanets. Following Bouc hv et al.l 
| |201Q) that found no evidence for a significant orbital eccentricity 
of WASP-21b we adopt a circular orbit. We included the quadratic 
limb dark ening (LD) coef ficients for the RISE filter V-l-R from the 
models of lHowarthl ( l20TTI) : a = 0.45451 and b = 0.210172. These 
were calculated for Te// = 5800 K log g=4.2 and [M/H ] =-0.5 
to match the stellar parameters from lSouchv et al] ( l20I(]h . We ini- 
tially kept the limb darkening parameters fixed during the fit. For 
each light curve, we included two extra parameters to account for a 
linear normalization. Therefore, 12 parameters were fitted. Besides 
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the linear normalization, no extra trends were removed from the 
light curve. 

To obtain the best fit para meters and uncertaint ies, we used a 
MCMC algorithm (e.g. | Tegmark et al. 2004: Collier Cameron et al.l 
l2007l : lGibson et al.ll2008l) . We begin by calculating the statistic 
of a set of proposed parameters, 



Table 1. Time residuals from the linear ephemeris. 



JV 

E 



(1) 



where fj is the flux observed at time j, rrij is the model flux and 
aj is the uncertainty of each fj as described in Section [J!2l At 
each step in the MCMC chain, each proposed parameter is per- 
turbed by a random amount which we call a "jump function". Each 
jump function is proportional to the uncertainty of each parameter 
multiplied by a random Gaussian number with mean zero and unit 
standard deviation. The new parameter set is accepted with proba- 
bility. 



P — rain 1, exp 



(2) 



where Ax'^ is the difference in the of subsequent parameters 
sets. Note, the new parameter set is always accepted if its is 
lower than the previous parameter set {P = 1). The jump functions 
are scaled by a common factor in order to rasure thal^25% of the 
steps are accepted, as suggested by Tegmark et alj j2004) . To esti- 
mate the uncertainty of each parameter and calculate the jump func- 
tions, an initial MCMC fit was performed. With these jump func- 
tions, we computed seven MCMC chains each of 150 000 points 
and different initial parameters. The initial 20% of each chain that 
corresponded to the burn in phase were discarded and the remaining 
parts merged into a master chain. We estimated the best fit param- 
eter as the mode of its probability distribution and the 1 a limits 
as the value at which the integral of the distributi on equals 0.341% 
from b oth sides of the mode. We computed the lOelman & Rubir] 
( Il992h statistic for each fitted parameter and concluded that chain 
convergence was good. 

To test how the limb darkening coefficients affect the derived 
system parameters, we repeated the MCMC procedure also fitting 
for the linear LD "a" which is the most sensitive to the observ- 
ing filter. The quadratic LD coefficient, h = 0.2 10172, was kept 
fixed, because as reported bv toibson et al.l j200 8^. the high preci- 
sion of the RISE light curves is not enough to fully constrain the 
LD coefficients (i.e., the MCMC does not converge when fitting 
both coefficients). We restricted the linear LD coefficient to be the 
same for all the light curves since they were all taken with the same 
filter. Therefore, in the second MCMC procedure we fitted 13 pa- 
rameters. We estimated a — 0.337 ± 0.034. 



4 RESULTS 

Comparing the fitted and fixed LD solutions, we concluded that 
although the fitted linear LD coefficient is statistically significantly 
different from the theoretical value, this does not affect the derived 
system parameters to any extent. The two solutions are within 1.5a 
of each ot her. Contrary to what was found by other authors (e.g., 
iGibson e t al. 2008; Southworth 2008), the uncertainties of the fitted 
LD solution are slightly smaller than those of the fixed solution. 
The of the fitted LD solution is similar to the fixed LD solution 
which does not justify the addition of an extra free parameter in the 
fit. Consequently, we conclude that our light curve is of insufficient 



Epoch Time residuals (sec) Uncertainty (sec) 



-78 


102 



11 
-7 



40 
24 
30 



Table 2. WASP-21 system parameters derived from the mcmc 



Parameter 



Value 



Normalised sepai'ation a/ij* 
Planet/star radius ratio Rp/Rt 
Orbital inclination / [degrees] 
Impact parameter b[R„] 
Transit duration Tt [days] 
Stellar density [p©] 



9 68+"-^'' 

87.34 ±0.29 
458+°-°''^ 
1430+0 0013 

652+0-041 



quality to better constrain the linear LD relative to that achieved by 
theoretical models and we choose to present the fixed LD solution. 

The estimate d transit times, combined the original ephemeris 
( iBouchv et aLll2O10i) were used to update the linear ephemeris. 



Tt{HJD) = T{0) + EP. 



(3) 



We found P = 4.3225060 ± 0.0000031 and To = 
2455084.51974 ± 0.00020 which was set to the mid transit time of 
the 2009 light curve. This ephemeris was used in the final MCMC 
procedures. 

For future reference, the time residuals from the linear 
ephemeris are given in Table [T] We conclude that the time resid- 
uals of WASP-21b are consistent with a linear ephemeris. 

The geometric system parameters of WASP-21 and the la un- 
certainties derived from the MCMC analysis with fixed limb dark- 
ening coefficients are given in Table |2] These parameters are di- 
rectly measured from the transit light curve and are only weakly 
dependent of stellar properties through the limb darkening co- 
efficients. Note, all the derived parameters presented in Table |2] 
were calculated at each point of the chain. Therefore, the final 
derived values and errors were determined from their probability 
distribution as done for the fitted values. We obtain a significantly 
lower density than was previous reported in the discovery paper 
p, = 0.84±0.09pq). 



4.1 Stellar mass and age 

To obtain the stellar and planetary physical properties, the geomet- 
ric parameters have to be scaled with the stellar mass. The new 
high quality transit light curves give a direct estimate of the stel- 
lar density. This allows a more accurate estimation of stellar mass 
than log g derived from spectral analysis ( Sozzetti et al. 2007). Cur- 
rently there are two main methods to derive the stellar mass from 
the stellar dens ity. The first uses is ochrones and mass tracks from 
stellar models jSozzetti et al ]|2007h and the second us es an empir- 
ical calibration derived from stellar eclipsing binaries jTorres et al.l 
^010; Enoch et al. 2010). 

iBouchvetalJ ( I2OI0I) derived the stellar ma ss through the em- 
pirical calibration between Tsff,pt and [Fe/H ] dTorres et al.l20I(]|) 
with the parametrisation of lEnoch et alj | [201oh . Following the same 
procedure, with the improved p,, we derive a stellar mass of 




Figure 3. Isochrone models (solid lines) from iDemarg ue et al.' i2004j) 
for WA SP-21 using [Fe/H] =-0.47 and [M/H] =-0.4 from Bo uchv et 

The age in Gyr is marked in the left of the respective model. We also 
show the mass tracks (dashed lines) for a stellar mass of 1.0, 0.9, 0.86 and 
a8M0. We overplot the T^j j = 5800 value adapted from lBouchv et alj 
i2010l) and the new (p*/p0)-l/3. 



1.02 ± O.O5M0. In Table [3] we present the mass and radius of 
WASP-21 and WASP-21b and the la uncertainties derived from 
the MCMC for a stellar mass of l.O2±O.O5M0 . We obtain a signif- 
icantly larger stellar and planetary radius than previously reported. 
This is due to the main-sequence assumption in the previous anal- 
ysis which as discussed below is found to be invalid. 

We also estimate the stellar mass from stellar models by inter- 
polatin g the Yonsei- Yale stellar e volution tracks Demarque et al.l 
(I2OO4 ') using the metallicity from Bouchv et al] ( I2OI0 I). These evo- 
lution tracks are plotted in Figure [5] along with the position of 
WASP-21. From the isochrones, we estimate a lower mass of 
O.86±O.O4M0 and an age of 12±2 Gyr for WASP-21. In Figureg] 
we also show the evolutionary tracks for stellar masses of 1.0, 0.95, 
0.86 and 0.8 M0 adapted from lOemargue et alj [2004 ). These sug- 
gest that WASP-21 is close to, or is already in the hydrogen-shell 
burning phase and hence is evolving off the main-sequence. This 
implies that the assumption of a main sequence mass-radius rela- 
tionship in the original analysis of lBouchv et alj feOlOh is faulty. 

There is a significant difference between the mass derived 
from evolutionary models, A/* — O.86±O.O4M0 , and the mass de- 
rived fro m the empirical cali bration, M* — 1.02 ± O.O5M0 . In the 
past, the [Torres et al.l j201Ct) calibration was found to be in agree- 
ment and a more straight-forward alt ernative to the stellar models 
jTorres etalj|20ld : lEnoch et alj|2010 ) . Moreover, it has the advan- 
tage that it can be directly included in a transit fitting procedure 
tenoch et alj|201ol) . However, recently the same discrepancy be- 
tween empirical and isochrone masses was also found for WASP- 
37 (S impson et al. 2011) and WASP-39 (Faedietal. 2011). The 
[Torres et al.l ((2010) eclipsing binaries sample used for calibrating 
their relationship does not contain many low-metallicity systems, in 
particular in the low -mass regime. Therefore, this suggests that the 
iTorres et all ( l2010h calibration might not hold for metal poor stars 
specially in the low mass regime. For these reasons, for WASP-21, 
we favour the lower mass derived from the evolution models. For 
a stellar mass of M, — 0.86 ± O.O4M0 we present the stellar and 
planetary radii for the WASP-21 system in Table [3] along with the 
la uncertainties derived from the MCMC analysis. 



Figure 4. Evolutionary mass tracks from iDemarque et alj (2004') for the 
same stellar parameters as Fig.|3]from left to right for stellar masses of 1.0, 
0.95, 0.86 and 0.8 M0 . WASP-21 position is in the 0.86 Mq evolutionary 
mass track for which we also show the 8, 10, 12, and 14 Gyr points. 



Table 3. WASP-21 system stell ar and planetary par ameters derived using ei- 
ther th e empirical caUb ration of lTorres et al.l j2010l) or the YY stellar inodels 
iDemaraue et alj2004l) . 



Torres inodels 



YY models 



Stellar mass M* [Mq] 
Stellar radius _R, _R, [R0] 
Stellar surface gravity log [cgs] 
Orbital semimajor axis a [AU] 

Planet mass Mp [Mjup] 
Planet radius Rp [Rjup] 
Planet density pp [p j] 
Planet surface gravity log gp [cgs] 



1.02 ±0.05 

1 161+°°^'' 

1.101_g 

4.32 ±0.02 
0.052 ±0.001 

0.30 ±0.01 

1 9if,+0.048 
„ +0.014 

'^•^'^-o.ois 
2.71 ±0.02 



0.86 ± 0.04 

4.29 ±0.02 
0.0494 ± 0.0009 



0.27 ±0.01 
1 140+0.045 

n 181+0015 
U-i»i_o.020 

2.71 ± 0.02 



To summarise, we derive a lower stellar mass, O.86±O.O4M0 , 
and a lower planetary mass, 0.27 ± O.OlMjup. We estimate the 
inclination of the orbit to be 87.3 ± 0.3 degrees. The radius of 
the star is found to be 1.10 ± 0.03 R0 and the planet radius is 
1.14±0.04 Rjup ,yieldingaplanetarydensity of 0.18 ±0.02 Pj . 



4.2 Eccentricity 

iBouchv et all ilOld) found the eccentricity to be statistically in- 
distinguishable from zero, i.e. the does not significantly im- 
prove when adding the two additional parameters to the circular 
model. For these cases, allowin g the eccentricity to flo at tends 
to overestimate the eccentricity dLucv &Sweenevlll97lh . Hence, 
Bouchv et al. (2010) adopted a circular orbit. Assuming a tidal 
dissipation parameter between 10^ and 10®, the circularisation 
timescale for WASP-21b is approximately between 0.017 and 
0.17 Gyr, respectively. Since this is much shorter than the derived 
age for the system we expect a circular orbit. However, if the or- 
bit not circular assuming a zero eccentricity results in underesti- 
mated uncertainties. Therefore, it is interesting to investigate the 
effect of a small non-zero eccentricity in the system parameters 
and their uncertainties. As an example, we assume an eccentricity 
of 0.04 ± 0.04 which is consistent with the discovery paper. We 
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repeated the MCMC procedure allowing the eccentricity to float. 
Because the transit light curves do not constrain the eccentricity, 
we include a prior on the eccentricity of the form: 

(ecc - eccoY 

where we assume ecco = a^cc ~ 0.04. This prior is added to 
equation 1 at each step of the chain. From the posterior eccentricity 
distribution we obtain an eccentricity of 0.038 ± 0.036 which is 
close to the input value. 

The maximum effect of the eccentricity upon the derived pa- 
rameters corresponds to the case where the transit occurs close to 
periastron (cj = 90°) or apastron (cj — —90°). Hence, in order 
to investigate the maximum deviation from a circular orbit we as- 
sume uj — 90° . For this particular case by assuming a circular orbit 
we would be overestimating a/R,, inc and p*, and underestimat- 
ing the stellar and planetary masses and radii. The opposite would 
have happen if we have assumed uj = —90°. 

The derived eccentric solution is within one sigma of the cir- 
cular solution and the uncertainties of a/i?*, / and p* are ~ 30% 
larger. This results in an increased uncertainty of ~ 30% on the 
radii and ^ 20% on the masses. Hence, we conclude that if the or- 
bit eccentricity is < 0.038 the system parameters would be within 
~ 1.3cr of the values given in Tables 2 and 3. 



5 DISCUSSION AND CONCLUSION 

We have presented two high quality transit light curves of WASP- 
21b taken with RISE. Together with the previous RISE partial tran- 
sit, these were fitted with an MCMC procedure to update the pa- 
rameters of the system. We have been conservative in our error es- 
timates by scaling the and by including time correlated noise in 
our analysis. 

The derived stellar density p* = 0.65 ± 0.05po and the esti- 
mated age for the system, 12 ± 2 Gyr, suggest that WASP-21 is in 
the process of evolving off the main sequence. Therefore, the main- 
sequence mass-radius relation assumed for WASP-21 in the discov- 
ery paper was invalid which led to a significant overestimation of 
the stellar density, thus affe cting the derived plane tary properties. 
Using the stellar models o f iDemarque et alj j2004l) . we derived a 
significantly lower stellar, M* = 0.86 ± O.O3M0, and planetary 
mass, Mp — 0.27 ± O.OlMjup. This lower host star mass some- 
what compensates the lower stellar density which re sults in a stel- 
lar rad ius which is within la of the one presented bv lBouchv et al.l 
( l2010h . 

We obtained a slightly larger planetary radius, Rj, = 1.14 ± 
0.04 R jup, for WASP-21b than previously reported. iFortnev et al.l 
( |2007|) hydrogen/helium coreless models predict a radius of ~ 
1.06 R© which is consistent within 2a with our estimated ra- 
diu s without the n e ed for any extra heating mechanism. Follow- 
ing iLaughlin et al.l ( l201ll) we compute a radius anomaly, = 
0.09, for WASP-21b . This supports the correlation reported by 
iLaughlin et"d] j201lh . i.e. = T^^t, where Tequ is the equilib- 
rium temperature o f the planet, which is ~ 1320 K for WASP- 
21b. iBouchv et all ilOVj) argued that the density of WASP-21b 
strengthens the correlation between planetary density and host star 
metallicity for hot Saturns (Guillot et al. 2006). With the addi- 
tion of the latest Saturn-mass planet disco veries (e.g. WASP-39, 
lFaedietal.ll201lL WASP-40. lAnderson et al]|201 ih this correlation 
appears weaker. However, if we scale for the equilibrium temper- 
ature with, for example J?, the correlation with metallicity is still 



strong (see Figure 6 in lFaedietal.ll20TTh . Moreover, the corre- 
lation also holds fo r the more massive planets (see Figure 3 in 
ILaughlin etalj201lh . 

Exoplanet transit light curves are often affected by system- 
atic noise that can in some cases dominate the photometric noise. 
Therefore, it is important to minimise the sources of systematic 
noise. In Section 3.1, we show an example of systematic noise 
present in our exoplanet transit observations and suggest that the 
first step to decrease this noise is to maintain the star in the same 
pixel position in the CCD during the observations. We confirm that 
defocused observations can also help decreasing systematic noise, 
as well d ecreasing deadtime loss es and hence improving the signal- 
to-noise jSouthworth et alj |2009). The systematic noise in our ob- 
servations was due to the variation of the stellar position across the 
CCD. 
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